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Project Goals 

The objective of this research task is to experimentally and computationally study the unsteady 
dynamics of a liquid hydrogen wave front translating down a pipeline during the chill down mode and 
develop a comprehensive computational model to predict the associated flow fields, thermal fields, and 
residence time. In order to accomplish the objective, the following tasks are planned: a) experimentally 
study the dynamics of a propagating liquid hydrogen wave front, b) experimentally measure the heat 
transfer rate associated with film boiling beneath the wave front, c) experimentally investigate transient 
nucleate boiling heat transfer associated with liquid film flow, d) develop a comprehensive computational 
model to predict the flow and temperature fields associated with propagating liquid hydrogen waves, and 
e) develop engineering models for cryogenic practitioners to predict the required chill down residence 
time under a variety of different operating conditions. 

Accomplishments 

Abstract 

A pseudo-steady model has been developed to predict the chilldown history of pipe wall temperature 
in the horizontal transport pipeline for cryogenic fluids. A new film boiling heat transfer model is 
developed by incorporating the stratified flow structure for cryogenic chilldown. A modified nucleate 
boiling heat transfer correlation for cryogenic chilldown process inside a horizontal pipe is proposed. The 
efficacy of the correlations is assessed by comparing the model predictions with measured values of wall 
temperature in several azimuthal positions in a well controlled experiment by Chung et al. (2004). The 
computed pipe wall temperature histories match well with the measured results. The present model 
captures important features of thermal interaction between the pipe wall and the cryogenic fluid, provides 
a simple and robust platform for predicting pipe wall chilldown history in long horizontal pipe at 
relatively low computational cost, and builds a foundation to incorporate the two-phase hydrodynamic 
interaction in the chilldown process. 


Introduction 

The cryogenic chilldown is encountered in many applications but is of particular importance in 
cryogenic transportation pipelines. For example, in rockets or space shuttle launch facility, cryogenic 
liquids as fuel are filled from the storage tank to the internal fuel tanks of a space vehicle through a 
complex pipeline system. To avoid evaporated fuel entering the space vehicle, a cryogenic chilldown 
prior to the filling is required to reduce the pipe wall temperature to the saturation temperature of the 
cryogenic liquid. 

Cryogenic chilldown involves complicated hydrodynamic and thermal interactions among liquid, 
vapor, and solid pipe wall. There exist few basic experimental studies and modeling efforts for chilldown 
of cryogenic fluids. Studies on cryogenic chilldown started in 1960’s accompanying the development of 
rocket launching system. Early experimental studies were conducted by Burke et al. (ref. 1), Graham 
(ref. 2), Bronson et al. (ref. 3), Chi and Vetere (ref. 4), Steward (ref. 5) among other researches. Bronson 
et al. (ref. 3) studied the flow regime in a horizontal pipe during the chilldown by liquid hydrogen. The 
results revealed that the stratified flow is prevalent in the cryogenic chilldown. 
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Flow regimes and heat transfer regimes in the horizontal pipe chilldown were also studied by Chi and 
Vetere (ref. 4). Information of flow regimes was deduced by studying the fluid temperature and volume 
fraction during the chilldown. Several flow regimes were identified: single phase vapor, mist flow, slug 
flow, annular flow, bubbly flow, and single -phase liquid flow. Heat transfer regimes were identified as: 
single-phase vapor convection, film boiling, nucleate boiling, and single-phase liquid convection. 

Recently, Velat et al. (ref. 6) systematically studied cryogenic chilldown in a horizontal pipe. Their 
study included: visually recording the chilldown in a transparent Pyrex pipe, which is used to identify the 
flow regime and heat transfer regime; collecting temperature histories at different positions of wall in the 
chilldown; recording pressure drop along the pipe. Chung et al. (ref. 7) conducted a similar study on the 
nitrogen chilldown at relatively low mass flux and provided the data needed to assess various heat transfer 
coefficients in the present study. 

Burke et al. (ref. 1) developed a crude chilldown model based on the 1-D heat transfer through the 
pipe wall and the assumption of infinite heat transfer rate from the cryogenic fluid to the pipe wall. The 
effects of flow regimes on the heat transfer rate were neglected. Graham et al. (ref. 2) correlated heat 
transfer coefficient and pressure drop with the Martinelli number (ref. 8) based on their experimental data. 
Chi (ref. 9) developed a 1-D model for energy equations of the liquid and the wall, based on film boiling 
heat transfer between the wall and fluid. An empirical equation for predicting chilldown time and 
temperature was proposed. 

Steward (ref. 5) developed a homogeneous flow model for the chilldown. The model treated the 
cryogenic fluid as a homogeneous mixture. The continuity, momentum and energy equations of mixture 
were solved to obtain density, pressure and temperature of mixture. Various heat transfer regimes were 
considered: film boiling, nucleate boiling and single-phase convection heat transfer. Separate treatment of 
different heat transfer regimes resulted in a significant improvement in the prediction of the chilldown 
process. Homogeneous mixture model was also employed by Cross et al. (ref. 10) who obtained a 
correlation for the wall temperature in the chilldown with an oversimplified heat transfer model of the 
heat transfer between the wall and the fluid. 

Stratified flow regime, which is the prevalent flow regime in horizontal chilldown, was first studied 
by Chen and Banerjee (refs. 11 to 13). They developed a separated flow model for the simulating cool- 
down by a stratified flow in a hot horizontal pipe. Both phases were modeled using 1-D mass and 
momentum conservation equations. The wall temperature was computed using a 2-D transient heat 
conduction equation. Their prediction for wall temperature agreed well with their experimental results. 
Although a significant progress was made on handling the momentum equations, the heat transfer 
correlations employed were not as advanced. 



Figure 1 . — Schematic of chilldown heat transfer regimes in a 
horizontal pipe. 
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Typical chilldown involves several heat transfer regimes as shown in figure 1. Near the liquid front is 
the film boiling. The knowledge of the heat transfer in the film boiling regime is relatively limited, 
because i) film boiling has not been the central interest in industrial applications; and ii) high temperature 
difference causes the difficulties in experimental investigations. For the film boiling on vertical surfaces, 
early work was reported by Bromley (ref. 14), Dougall and Rohsenow (ref. 15) and Laverty and 
Rohsenow (ref. 16). Film boiling in a horizontal cylinder was first studied by Bromley (ref. 17); and the 
Bromley correlation was widely used. Breen and Westwater (ref. 18) modified Bromley’s equation to 
account for very small tubes and large tubes. If the tube is larger than wavelength associated with the 
Taylor instability, the heat transfer correlation is reduced to Berenson’s correlation (ref. 19) for a 
horizontal surface. 

Empirical correlations for film boiling were proposed by Flendrick et al. (refs. 20 and 21), Ellerbrock 
et al. (ref. 22), von Glahn (ref. 23), Giarratano and Smith (ref. 24). These correlations relate a simple or 
modified Nusselt number ratio to the Martinelli parameter. Giarratano and Smith (ref. 24) gave detailed 
assessment on these correlations. All these correlations are for steady state cryogenic film boiling. They 
may not be suited for transient chilldown application. 

When the pipe chills down further, film boiling ceases and transient boiling occurs, followed by 
nucleate boiling. The heat transfer in the transient boiling is even more complicated and it is usually 
assumed that the boiling switches from film boiling to nucleate boiling right away. The position of film 
boiling transitioning to the nucleate boiling is often called rewetting front, because from that position the 
cold liquid starts touching the pipe wall. Usually the Leidenfrost temperature indicates the transition from 
the film boiling to the nucleate boiling. 

The study on the convection nucleate boiling is extensive. A general correlation for saturated boiling 
was introduced by Chen (ref. 25). Gungor and Winterton (ref. 26) modified Chen’s correlation and extend 
it to subcooled boiling. Enhancement and suppression factors for macro-convective heat transfer were 
introduced. Gunger and Winterton’ s correlation can fit experimental data better than the modified Chen’s 
correlation (ref. 27) and Stephan and Auracher correlation (ref. 28). Kutateladze (ref. 29) and Steiner (ref. 
30) also provided correlations for cryogenic fluids in pool boiling and forced convection boiling. 

Although they are not widely used, they should be more applicable for cryogenic fluids as the correlation 
was directly based on cryogenic conditions. 

As the wall temperature drops further, the nucleate boiling is replaced by pure convection. The 
convection heat transfer can be modeled using Dittus-Boelter equation (ref. 31) for the fully developed 
turbulent pipe flow or corresponding laminar heat transfer equation for laminar pipe flow (ref. 31). The 
condition in which nucleate boiling switches to single-phase convection heat transfer is that all nucleate 
sites are suppressed. 

Although two-fluid model can describe the fluid dynamics aspect of the chilldown process, it suffers 
from computational instability for moderately values of slip velocity between two phases which limits its 
application. To gain fundamental insight into the thermal interaction between the wall and the cryogenic 
fluid and to be able to rapidly predict chilldown in a long pipe, an alternative pseudo-steady model is 
developed. In this model, liquid wave front speed is assumed to be constant and is the same as the bulk 
liquid speed (ref. 32). It is also assumed that steady state thermal fields for both the liquid and the solid 
exist in a reference frame that is moving with the wave front. The heat transfer between the fluid and the 
wall is modeled using different heat transfer correlations depending on the operating heat transfer regime 
at a given location. Various improvements on the correlations are introduced, including the development 
of a new film boiling heat transfer coefficient. The governing equation for the solid thermal field becomes 
a parabolic equation that can be efficiently solved. It must be emphasized that a great advantage of the 
pseudo-steady model is that one can assess the efficacy of the film boiling model independently from that 
of the nucleate boiling model since the down stream information in the nucleate boiling regime cannot 
affect the temperature in the film boiling regime. In another word, even if the nucleate boiling heat 
transfer coefficient is inadequate, the film boiling heat transfer coefficient can still be assessed in the film 
boiling regime by comparing with the measure temperature for the right period of time. After the 
satisfactory performance is achieved for the film boiling regime, the nucleate boiling heat transfer model 
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can be subsequently assessed. In the Results section, those detailed assessments of the heat transfer 
coefficients are provided by comparing the computed temperature variations with the experimental 
measurements of Chung et al. (ref. 7). Satisfactory results are obtained. 
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Subscripts 

0 

i and o 


Boiling number 
solid heat capacity 
heat capacity 
pipe diameter 
thickness of pipe wall 
gravity 

heat transfer coefficient 

pool boiling heat transfer coefficient 

convection boiling heat transfer coefficient 

latent heat 

Jacob number 

thermal conductivity 

effective thermal conductivity 

Nusselt number 

pressure 

radius of pipe 

inner and outer radius of pipe 
Rayleigh number 
Reynolds number 
Peclet number 
Prandtl number 
suppression factor 
temperature 

Leidenfrost temperature 

transition temperature between nucleate boiling to convection heat transfer 

room temperature 

time 

velocity 

vapor film velocity 
vapor film coordinates 
transformed coordinate 
cylindrical coordinates 
liquid volume fraction 
Martinelli number 
vapor film thickness 
emissivity 

azimuthal coordinate 
viscosity 

dimensionless temperature; azimuthal coordinate 
liquid surface tension; Stefan Boltzmann constant 


characteristic value 
inner and outer pipe 
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I liquid 

v vapor 

w wall 

sat saturated 

Superscripts 

dimensionless variable 


Formulation 

In pseudo-steady chilldown model, it is assumed that both the liquid and its wave front moves at a 
constant speed U which is taken from estimated experimental condition in the present model. Thus the 
main emphasis of the present study is on the modeling of the heat transfer coefficients in different heat 
transfer regimes and the computation of the thermal field within the solid pipe. 

Solid Heat Transfer 

The thermal field inside the solid wall is governed by the 3-D unsteady energy equation: 


pc 


dT = d_ 
dt dz 
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( 1 ) 


Since the wave front speed U is assumed to be a constant, it can be expected that when the front is 
reasonably far from the entrance region of the pipe, the thermal field in the solid is in a steady state when 
it is viewed in the reference frame that moves with the wave front. Thus, the following coordinate 
transformation is introduced, 


Z = z + Ut . 


( 2 ) 


Equation (1) is then transformed to: 
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For further simplification, the following dimensionless parameters are introduced, 


( 3 ) 



Figure 2. — Coordinate systems: laboratory frame is 
denoted by z; moving frame is denoted by Z. 
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where k 0 is the characteristic thermal conductivity, and Co is the characteristic heat capacity. Equation (3) 
can be normalized as 
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where Pc = 

ko 


is the Peclet number. 


Under typical operating condition for the cryogenic chilldown process, Pc ~ 0( 1 0 2 ) - 0(19 3 ). The 
first term on the RHS of eq. (5) is small compared with the rest of terms and thus can be neglected, eq. (5) 
becomes 
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( 6 ) 


which is parabolic. Hence in the Z'-direction, only one condition is needed. In the (p-direction, periodic 
boundary conditions are used. On inner and outer surfaces of the pipe wall, where r = R] and R >, proper 
boundary conditions for the temperature are required. 

For convenience, Z' = 9 is set at the liquid wave front. In the region of Z! < 9, the inner wall is 
exposed to the pure vapor. Although there may be some liquid droplets in the vapor that can cause 
evaporative cooling when the droplets deposit on the wall and the cold flowing vapor can absorb some 
heat from the wall, the heat flux due to these two mechanisms is small compared with the heat transfer 
between liquid and solid wall in the region of Z' > 9. Hence, heat transfer for Z' < 9 is neglected and it is 
assumed that 9 = 1 at Z' = 9. The computation starts from the Z! = 9 to Z'— until a steady state solution 
in the Z'-direction is reached. An implicit scheme in the Z'-direction is employed to solve eq. (6). 

Liquid and Vapor Flow 

The two-phase flow is assumed to be stratified as was observed in (ref. 7). Both liquid and vapor 
phases are assumed to be at the saturated state. Liquid volume fraction determining which part of wall is 
in contact with the liquid or the vapor, is specified at every cross-section along the Z'-direction based on 
experimental information. The visual studies (refs. 6 and 7) show that the liquid volume fraction increases 
gradually, rather than abruptly, near the liquid wave front and becomes almost constant during most of the 
chilldown. Hence, the following liquid volume fraction a as a function of time is used in the computation 
of the solid-fluid heat transfer coefficient, 


/ 

a = oco s i n 
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t 

f 0 


t<t 0 


a = a 0 


t>t 0 


(7) 


where to is characteristic chilldown time, and oco is characteristic liquid volume fraction. Here the time 
when the nucleate boiling is almost suppressed and slope of the wall temperature profile becomes flat is 
set as characteristic chilldown time. 

The vapor phase velocity is assumed to be a constant. However, it was not directly measured. It is 
computationally determined by trial-and-error by fitting the computed and measured wall temperature 
variations for numerous positions. 
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Heat Transfer Between Cryogenic Fluid and Solid Pipe Wall 

During the chilldown, the fluid in contact with the pipe wall is either liquid or vapor. The mechanisms 
of heat transfer between the liquid and the wall and between the vapor and the wall are different. Based 
on experimental measurements and theoretical analysis, liquid-solid heat transfer accounts for a majority 
of the total heat transfer. However, the prediction for heat transfer is much more complicated than the 
heat transfer between the vapor and the wall. The heat transfer between the liquid and the wall is 
discussed first. 

Heat Transfer Between Liquid and Solid wall 

Heat transfer between liquid and solid wall includes film boiling, nucleate boiling, and single-phase 
convection heat transfer. The transition from one type of boiling to another depends on many parameters, 
such as wall temperature, wall heat flux, and various properties of fluid. For simplicity, a fixed 
temperature approach is adopted to determine the transition points. That is, if the wall temperature is 
higher than the Leidenfrost temperature, film boiling dominates. If the wall temperature is between 
Leidenfrost temperature and a transition temperature, T 2 , nucleate boiling takes place. Here, the transient 
boiling regime between film boiling and nucleate boiling is neglected. The reason is mainly the difficulty 
associated with the determination of the conditions for the transient boiling to occur and the lack of a 
reliable transient heat transfer correlation. If the wall temperature is below the transition temperature T 2 , 
pure convection heat transfer dominates. The values of Leidenfrost temperature and transition 
temperature are determined by matching model prediction with the experimental results. 

Film boiling heat transfer . — Because of the nature of high wall superheat in chilldown, the film 
boiling plays a major role in chilldown in terms of the time span and in terms of the total amount of heat 
removed from the wall. Currently there exists no specific film boiling correlation for chilldown 
applications with such high superheat. If one uses conventional film boiling correlations, necessary 
modifications for cryogenic application must be made for chilldown. 

One of cryogenic film boiling heat transfer correlations was provided by Giarratano and Smith 
(ref. 24) 


Nu 

. Nu ca i c 


*Bo -°' 4 =/(X«) 


( 8 ) 


where Nu ca / c is the Nusselt number for the forced convection heat transfer. In this correlation, the heat 
transfer coefficient is the averaged value for the whole cross section. Similar correlations for cryogenic 
film boiling also exist in literature. The correlations were obtained from measurements conducted under 
steady states. The problem with the use of this steady state film boiling correlation is that it does not take 
into account the change of flow regime as encountered in the chilldown. For example, for the same 
quality, the heat transfer rate in annular flow is much different from that in stratified flow, while those 
empirical correlations cannot take such difference into account. 

Furthermore, in this study, local heat transfer coefficient is needed in order to incorporate the thermal 
interaction with the pipe wall. Since the two-phase flow regime information is available through 
visualization in the present study, the modeling effort needs to take into account the knowledge of the 
flow regimes. 

There are several correlations for the film boiling based on the analysis of vapor film boundary layer 
and stability of the thin vapor film, such as Bromley’s correlation (ref. 17) and Breen and Westerwater’s 
correlation (ref. 18) for film boiling on the outer surface of a hot tube, Frederking and Clark’s (ref. 33) 
and Carey’s (ref. 34) correlations for film boiling on the surface of a sphere. However, none of these was 
obtained for cryogenic fluids or for film boiling on the inner surface of a pipe. 

A new correlation for the film boiling in cryogenic chilldown inside a tube is presented here. The 
schematic of film boiling inside a pipe is shown in figure 3 with a cross-sectional view. Bulk liquid is 
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vapor 



Figure 3. — Schematic diagram of film boiling in stratified 
flow. 


near the bottom of the pipe. Beneath the liquid is a thin vapor film. Due to buoyancy force, the vapor in 
the film flows upward along the azimuthal direction. Heat is transferred through the thin vapor film from 
the solid to the liquid. Reliable heat transfer correlation for film boiling in pipes or tubes requires the 
knowledge of the thin vapor film thickness which can be obtained by solving the film layer continuity, 
momentum, and energy equations. 

To simplify the analysis for the vapor film heat transfer, it is assumed the liquid velocity in the 
azimuthal direction is zero and the vapor flow in the direction perpendicular to the cross-section is 
negligible. It is further assumed that the vapor film thickness is small compared with the pipe radius and 
vapor flow is in steady state, incompressible and laminar. The laminar flow assumption can be confirmed 
post priori as the Reynolds number, Re, based on the film velocity and film thickness is typically of 
<9(10° ~ 10°). In terms of the x- and v-coordi nates and (, u , v) velocity components shown in figure 3, the 
governing equations for the vapor flow are similar to boundary-layer equations: 
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where subscribe v refers to properties of vapor. 

Because the length scale in the azimuthal (x) direction is much larger than the length scale at the 
normal (y) direction, the v-component may be neglected. Furthermore, the convection term is assumed 
small and is neglected. The resulting momentum equation is simplified to 


1 dp 
p v dx 
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The vapor pressure can be evaluated by considering the hydrostatic pressure from liquid core as: 
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(13) 


p = Po + Pi gR 
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where (|) is angular position where the film merges with the vapor core. The momentum equation becomes 
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The vapor velocity boundary condition is u = 0 at y = 0 and u - ui - 0 at y = 8. The vapor velocity profile 
is: 
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The energy and mass balances on the vapor film requires that 
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Neglecting the convection, the vapor energy equation is: 
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dy 2 
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The following linear temperature profile is thus obtained, 
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Substituting the temperature and velocity profiles into eq. (17) yields 
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Equation (20) has analytical solution: 
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where Ja is Jacob number and Ra is Raleigh number: 
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and F(Q) 
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describes the geometric dependence of vapor film thickness. 
The mean velocity u as a function of 0 is thus 


(24) 
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Curves for F(Q) and F 2 (0)sin 0 based on numerical integration are shown in figure 4. The vapor film 
thickness has a minimum at 0 = 0 and is nearly constant for 0 < tt/2. It rapidly grows after 0 > n/2. The 
singularity at the top of tube when 0 — > 71 is of no practical significance since the film will merge with the 
vapor core at the vapor-liquid interface. The vapor velocity is controlled by F(0) 2 sin 0 which is zero at the 
bottom of the pipe and increases almost linearly in the lower part of the tube where the vapor film 



Figure 4. — Numerical solutions of the vapor thickness and velocity 
influence functions. 
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thickness does not change substantially. In the upper part of the tube, due to the increase in the vapor film 
thickness, the vapor velocity gradually drops back to zero at the top of the tube. Thus a maximum velocity 
may exist in the upper part of the tube. 

The local film boiling heat transfer coefficient is easily obtained from the linear temperature profile. 

It is 


h = ^- = 0.6389 
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k v f Ra 'i 4 

DF(d){ Ja J 


(25) 


Nucleate Boiling and Convection Heat Transfer 

For the nucleate flow boiling heat transfer, Gungor and Winterton’s correlation (ref. 26) is widely 
used due to its relatively better accuracy in predicting heat transfer coefficients. However, a closer 
examination on this correlation shows that it is based mainly on the following parameters: Pr, Re, and 
quality x. Similar to the development of conventional film boiling correlations, these parameters all reflect 
overall properties of the flow in the pipe and are not directly related to flow regimes. Thus it cannot be 
used to predict the local heat transfer coefficients in chilldown. 

Chen’s (ref. 25) flow boiling correlation is based on separating the heat transfer to micro- and macro- 
convection heat transfer. Micro-convection heat transfer represents the contribution from boiling heat 
transfer, and macro-convection represents the contribution from the forced convection heat transfer. 
However his correlation fits best for annular flow because of the derivation of suppression factor. In 
stratified flow regime, which is common in cryogenic chilldown, Chen’s correlation may not be 
applicable. 

Several correlations have been tried in this study, including Gungor and Winterton’s correlations 
(ref. 26), Chen’s correlation (ref. 25), and Kutateladze’s correlations (ref. 29). None gives a satisfactory 
heat transfer rate that is need to match the experimentally measured temperature histories in (ref. 7) in the 
nucleate boiling regime. Among them, Kutateladze’s correlation gives more reasonable result. In this 
correlation, the total heat transfer coefficient h is 


h h COm , + hpggl , 


where h com is given by Dittus-Boelter equation for fully developed pipe flow, 

h conv = 0.023* Re?' 8 Pr®' 4 * k t /D 
and pool nucleated boiling heat transfer coefficient h poo \ is 
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(27) 


(28) 


in which AT is wall superheat. 

Kutateladze’s correlation (ref. 29) was proposed without considering the effect of nucleate site 
suppression. This obviously leads to an overestimation of the nucleate boiling heat transfer rate. Hence a 
modified version of Kutateladze’s correlations (28) is used, 

h = h CO nv+S*h p gg,, (29) 

with S being the suppression factor and h poo i is given by (28). 
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When AT drops to a certain range all the nucleate sites are suppressed. The heat transfer is dominated 
by single phase forced convection. The heat transfer coefficient can then be predicated using Dittus- 
Boelter equation, eq. (27), if flow is turbulent, or eq. (30), if flow is laminar. 

h conv = 4.36* k,/D. (30) 

Heat Transfer between Vapor and Solid Wall 

The heat transfer between the vapor and the wall can be estimated by treating the flow as a fully 
developed convection flow, neglecting the liquid droplets that are entrapped in the vapor. The heat 
transfer coefficient of vapor forced convection flow is 

h,, = 0.023 * Re °- 8 Pr v 0 - 4 *k v ID (turbulent flow) (3 1 ) 

or 

h v =436*k v /D (laminar flow) (32) 

Heat Transfer between Solid Wall and Environment 

For a cryogenic flow facility, although serious insulation is applied, the heat leakage to environment 
is still considerable due to the large temperature difference between the cryogenic fluid and the 
environment. It is necessary to evaluate the heat leakage from the inner pipe to environment in order to 
make realistic assessment of the model prediction with the experimental results (ref. 7). 

Vacuum insulation chamber between the inner and outer pipes is used in cryogenic transport pipe 
(ref. 7), as shown in figure 5. Radiation heat transfer exists between the inner and outer pipe. 

Furthermore, the space between the inner and outer pipe is not an absolute vacuum. There is residual air 
that causes free convection between the inner and outer pipes driven by the temperature difference of the 
inner and outer pipe. 

The radiation between the inner pipe and outer pipe becomes significant when the inner pipe is 
chilled down. The heat transfer coefficient is proportional to the difference of the fourth power of wall 
temperatures. Exact evaluation of the radiation heat transfer between the inner and outer pipe is a difficult 
task. Hence a simplified model based on the overall radiation heat transfer between long concentric 
cylinders with constant temperature at inner pipe and outer pipe (ref. 31) is used to evaluate the heat 
transfer rate at every axial location of the pipe. It is not quantitatively correct, but can provide reasonable 
estimate for the magnitude of the radiation heat transfer between pipes through the vacuum. The local 
radiation heat transfer rate per unit area on the surface of inner pipe q' ra( j is 



Figure 5. — Schematic of vacuum insulation chamber. 
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where the a is Stefan Boltzmann constant, ( r„ 8,) and (r„, e 0 ) are the radius and emissivity of inner pipe 
and outer pipe, respectively, and T waU is the local inner wall temperature, T 0 is the room temperature that 
is assumed constant in the entire outer pipe. Here the emissivity is also assumed to be constant during the 
entire chilldown process. 

For the free convection heat transfer due to the residual air in the vacuum chamber between the inner 
pipe and outer pipe, Raithby and Hollands’ correlation (ref. 35) is used for the heat transfer rate. The 
average heat transfer rate per unit length of the cylinder is 




(34) 


where T is assumed constant on the inner and outer wall along the azimuthal directions, and k e jj is the 
effective thermal conductivity. 

Results and Discussions 

In experiment by Chung et al. (ref. 7), liquid nitrogen was used. The flow regime is revealed to be 
stratified flow by visual observation, as shown in figure 7, and wall temperature history in several 
azimuthal positions is measured by thermal couples and recorded on a computer. 

Experiment of Chung et al. 

In experiment of Chung et al. (ref. 7), a concentric pipe test section (fig. 6) was used. The chamber 
between the inner and outer pipe is vacuumed, but about 20 percent air remained. The inner diameter 
(i.d.) and outer diameter (o.d.) of inner pipe are 11.1 and 15.9 mm, and i.d. and o.d. of outer pipe are 95.3 
and 101.6 mm, respectively. Numerous thermal couples were placed at different locations of the inner 
pipe. Some were embedded very closely to the inner surface of the inner pipe while others measure the 
outside wall temperature of inner pipe. Experiments were carried out at room temperature and atmosphere 
pressure. Liquid nitrogen flows from a reservoir to the test section driven by gravity. As liquid nitrogen 
flows through the pipe, it evaporates and chills the pipe. Some of the typical visual results are shown in 
figure 7. The measured average liquid nitrogen velocity is U ~ 5 cm/s. Vapor velocity is not measured in 
this experiment. In this study, it is determined through trial-and-error by fitting the computed and 
measured temperature histories. The characteristic liquid volume fraction is 0.3 from the recorded video 
images. The characteristic time used in this computation is t 0 = 100 s. The Leidenfrost temperature for the 
nitrogen is around 1 80 K; hence the temperature when the film boiling ends and nucleates boiling starts is 
set at 1 80 K. The transition temperature at which the nucleate flow boiling switches to purely convection 
heat transfer is 140 K based on experimental results. The material of the inner pipe and outer pipe used in 
the experiment of (ref. 7) are Pyrex glass with emissivity of 0.82 at room temperature. 
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S’tondord to 



Figure 6. — Schematic of Chung et al.’s cryogenic two- 
phase flow test apparatus (ref. 7). 




wwrn 







Figure 7. — Experimental visual observation of Chung et al.’s 
cryogenic two-phase flow test (ref. 7). 


Comparison of Pipe Wall Temperature 

In the computation, there are 40 grids along the radial direction and 40 grids along azimuthal 
direction for the inner pipe (fig. 8). The results of temperature profile at 40x40 grids and higher grids 
resolution show that 40x40 grids are sufficient. Figures 9 and 10 compare the measured and computed 
wall temperature as a function of time at positions 11, 12, 14 and 15 shown in figure 8. For the modified 
Kutateladze’s correlation a suppression factor 0.01 is used. Vapor velocity is 0.5m/s based on the best fit. 
The overall temperature histories agree well in the film boiling stage. Thus, the film boiling heat transfer 
coefficient based on the first principle and incorporated the flow structure gives very reliable prediction. It 
must also be noted that the value of the Leidenfrost temperature does not affect the computed temperature 
history prior to that transition point since the governing equation is parabolic. The good agreement before 
the Leidenfrost temperature is reached is entirely due to the superior performance of the new film boiling 
heat transfer coefficient. 

During the stage of the rapidly decreasing wall temperature after the Leidenfrost temperature, the 
computed wall temperature drops slightly faster than the measured value. The rapid decrease in the wall 
temperature is due to initiation of nucleate boiling which is more efficient for heat transfer than the film 
boiling. Reasonable agreement between the computed and measured histories in this nucleate boiling 
regime is due to: i) the good agreement already achieved in the film boiling stage; ii) valid choice for the 
Leidenfrost temperature that switches the heat transfer regime correctly; and iii) appropriate modification 
of Kutateladze’s correlations. 
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In the final stage of chilldown, the wall temperature decreases slowly, and the computed wall 
temperature shows the same trend as the measured one but tends to be a little lower. Figure 1 1 shows the 
temperature distribution of a given cross-section at different times of chilldown. Because the upper part of 
pipe wall is exposed to nitrogen vapor, cooling effect is much reduced. 



Figure 8. — The computational grids arrangement and 
positions of thermal couples. 



Figure 9. — Comparison of wall temperature of 
transducer 12 and 15, which is at the bottom of pipe, 
between the experiment and computation. For 
comparison purpose, the solution in the Z-direction is 
converted to time t. 
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Figure 10. — Comparison of wall temperature of 
transducer 1 1 and 14, which are at the upper part of 
the pipe, between the experiment and computation. 
For comparison purpose, the solution in the Z- 
direction is converted to time t. 
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Figure 1 1 . — Wall temperature distribution across the section at t = 0, 50, 100 and 300 sec. 
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Discussions and Remarks 


In figure 9, the wall temperature based on the film boiling correlation of Giarratano and Smith 
(ref. 24) is also shown. Apparently, the correlation of Giarratano and Smith (ref. 24) gives a very low heat 
transfer rate so that the wall temperature remains high. This comparison confirms our earlier argument 
that correlations based on the overall flow parameter, such as quality and averaged Reynolds number, are 
not applicable for the simulation of unsteady chilldown. 

The nucleate flow boiling correlations of Gungor and Winterton (ref. 26), Chen (ref. 25), 
Kutateladze’s correlations (ref. 29) are also compared in this study. Gungor and Winterton’s correlation 
fails to give a converged heat transfer rate during the transition from the film boiling to the nucleate 
boiling. Chen’s correlation overestimates the heat transfer rate, and cause unrealistically large temperature 
drop on the wall, which results in the halt of the computation. Only Kutateladze’s correlation gives a 
reasonable heat transfer rate. However, the temperature drop near the bottom of the pipe is still faster than 
the measured one as shown in figure 9. This may be due to the fact that most of nucleate boiling 
correlations were obtained from experiments of low wall superheat. However, in this cryogenic 
chilldown, wall superheat is much higher than that in the normal nucleate boiling experiments. Another 
reason is that the original Kutateladze’s correlation does not include suppression factor. This leads to 
overestimating heat transfer coefficient. The modified Kutateladze’s correlation with suppression factor 
gives reasonable chilldown result in figure9. 

Further examination of figures 9 and 10 indicates that although we have considered the heat leak from 
the outer wall to the inner wall through radiation and free convection, the computed temperature is still 
lower than the measured temperature during the final stage of chilldown, where the heat transfer rate 
between fluid and wall is low due to lower wall superheat and possibly heat leak. The temperature 
difference between the computed and measured values at position 12 and position 15 suggests that there 
may be additional heat leak, which affects the measurements but is not taken in account in the present 
modeling. 

In this study, the pseudo-steady chilldown model is developed to predict the chilldown process in a 
horizontal pipe in stratified flow regime. This model can also be extended to describe annular flow 
chilldown in the horizontal or the vertical pipe with minor changes on the boundary condition for solid 
temperature. It can also be extended to study the chilldown in the slug flow as long as we specify the 
contact period between the solid and the liquid or the vapor. The disadvantage of the current pseudo- 
steady chilldown model is that the fluid interaction inside the pipe is largely neglected and both vapor and 
liquid velocities are assumed to be constant. Compared with a more complete model that incorporates the 
two-fluid model, the present pseudo-steady chilldown model requires more experimental measurements 
as inputs. However, the pseudo-steady chilldown model is computationally more robust and efficient for 
predicting chilldown process. It provides overall reasonable results for the solid wall temperature. While a 
more complete model for chilldown process that incorporates the mass, momentum, and energy equations 
of vapor and liquid is being developed to reduce the dependence of the experimental inputs for the liquid 
velocity and trial-and-error for the vapor velocity, the present study has revealed useful insights into the 
key elements of two-phase heat transfer encountered in the chilldown process which have been largely 
ignored. It also laid the necessary modeling foundation for the incorporation of the two-fluid model. 

Conclusions 

A pseudo-steady chilldown computational model has been developed to understand the heat transfer 
mechanisms of cryogenic chilldown and predict the chilldown wall temperature history in a horizontal 
pipeline. The model assumes the constant speed of the moving liquid wave front, and steady state of the 
thermal field in the solid in a reference frame that moves with the liquid wave front, and saturated 
temperature of both liquid and vapor so that the 3-D unsteady problem can be transformed to a 2-D, 
parabolic problem. A new film boiling heat transfer coefficient in the cryogenic chilldown condition is 
developed using the first principle and incorporating the stratified flow structure. The existing nucleate 
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boiling heat transfer correlations may not work well under cryogenic condition. A modified Kutateladze’s 
correlations with suppression factor adequately describes heat transfer coefficient. With the new and 
modified heat transfer correlations, the pipe wall temperature history based on the pseudo-steady 
chilldown model matches well with the experimental results by Chung et al. (ref. 7) for almost the entire 
chilldown process. The pseudo-steady chilldown model has captured the important features of thermal 
interaction between the pipe wall and the cryogenic fluid. 
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